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We show that the temperatures of the emergent non-electron neutrinos and the binding energy 
released by a galactic Type II supernova are determinable, assuming the Large Mixing Angle (LMA) 
solution is correct, from observations at the Sudbury Neutrino Observatory (SNO) and at Super- 
Kamiokande (SK). If the neutrino mass hierarchy is inverted, either a lower or upper bound can be 
placed on the neutrino mixing angle 0 13 , and the hierarchy can be deduced for adiabatic transitions. 
For the normal hierarchy, neither can 0 13 be constrained nor can the hierarchy be determined. Our 
conclusions are qualitatively unchanged for the proposed Hyper-Kamiokande detector. 


Neutrino oscillations convincingly explain the solar 
and atmospheric neutrino anomalies [1,2]. Atmospheric 
data and initial K2K data [3] indicate oscillations with 
IAto^ = m3 — m 2 | ~ 3 x 10 -3 eV 2 and sin 2 2023 ~ 1 [4] 
that will be tested to within 10% accuracy at MI¬ 
NOS [5,6]. (For the standard parameterization of the 
neutrino mixing matrix see Ref. [7]). The large mixing 
angle (LMA) solution (Am^ ~ 5 x 10~ 5 eV 2 , sin 2 2012 ~ 
0.8) which is emerging as the solution to the solar neu¬ 
trino problem [8] will be tested to 10% accuracy [9] at 
KamLAND [10]. Thus, in the near future, all parameters 
relevant to neutrino oscillations will be known, except 
sin 2 2013, which is bounded above by 0.1 at the 95% C.L. 
by the CHOOZ experiment [11], the sign of Am^ and 
the CP violating phase. In the longer term, long base¬ 
line neutrino experiments using upgraded conventional 
neutrino beams could achieve a sensitivity to sin 2 2013 
of about 10 -3 [12], but neutrinos from a galactic super¬ 
nova can probe values that are more than two orders of 
magnitude smaller. 

The objective of this Letter is to determine what 
the expected neutrino signals at SNO and SK from a 
Type II galactic supernova can tell us about sin 2 20i3, 
sgn(A?n 2 1 ), the neutrino temperatures, and the binding 
energy released in such an event if the LMA solution is 
confirmed. Throughout, we assume that solar and at¬ 
mospheric parameters will be known to within 10% from 
upcoming experiments. Since low energy and u T are 
indistinguishable at SNO and SK, only their transitions 
with u e can be studied. 

Supernova neutrinos: During the early stages of a su¬ 
pernova explosion, as the shock wave rebounds from the 
dense inner core of the star and crosses the electron neu- 


trinosphere, z/ e ’s from electron capture on protons are 
released resulting in a breakout or neutronization burst 
that carries away ~ 10 51 ergs. The duration of this 
burst lasts only a few milliseconds (no more than 10) 
and any non-electron neutrino events at SNO during this 
time are a consequence of v e —> z/ MiT oscillations. For 
progenitor stars of mass ~ 15M 0 , numerical simulations 
find that following the neutronization burst, 99% of the 
binding energy released, Ej, = 1.5 — 4.5 x 10 53 ergs, is 
roughly equipartitioned in the form of neutrinos and an¬ 
tineutrinos of all flavors [13]. Including effects of nu¬ 
cleon bremsstrahlung and electron neutrino pair anni¬ 
hilation, the luminosities are approximately related by 
L Ve ~ L De ~ (1-2 )L Vx where x = n,fb,r,f [14,15]. 
This emission occurs on a timescale of tens of seconds. 
The mean energies of the different flavors of neutri¬ 
nos are determined by the strength of their interactions 
with matter, with the most strongly interacting neutri¬ 
nos leaving the star with the lowest mean energy Le., 
(E Ue ) < (Ep e ) < (E Vx ). The authors of Ref. [14] (see also 
Ref. [16]) emphasize that spectral differences are very 
small, typically {E„ e ) : (E De ) : (E„ x ) :: 0.85 : 1 : 1.1. The 
spectra of neutrinos can be modeled by pinched Fermi- 
Dirac distributions. We can write the unoscillated differ¬ 
ential flux at a distance D from the supernova as 

F = _ ^ _^_ (l) 

“ 247rD 2 T4|Li 4 (-e^)| e E / T «~^ + 1 ’ 1 J 

where a = v e ,v e ,v x , Li n (z ) is the polylogarithm func¬ 
tion and r] a is the degeneracy parameter. The tem¬ 
perature of the neutrinos, T a , is related to (E a ) via 
(E a ) = 3 ^|~^ Q j r a . We shall use (E a ) and T a inter¬ 
changeably since they are equivalent to each other once 
rj a is specified. Strictly speaking, weak magnetism effects 
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may result in being about 7% higher than T v [17]. 
However, we have explicitly checked that the inequality 
of these temperatures does not affect our results. 

As the neutrinos leave the star, they encounter a den¬ 
sity profile that falls like 1/r 3 [18]. If the mass hierarchy 
is normal, i.e. Am§ 1 > 0, (inverted, i.e. < 0), 

neutrinos (antineutrinos) pass through a resonance at 
high densities (10 3 —10 4 g/cm 3 ) which is characterized by 
(Am 31 , sin 2 20i 3 ) and the neutrinos pass through a sec¬ 
ond resonance at low densities (~ 20 g/cm 3 for the LMA 
solution) that is determined by (Am| 1; sin 2 2 ^ 12 ) [19]. 
Transitions in the latter resonance are almost adiabatic, 
with an essentially zero probability of level crossing. We 
denote the jumping probability in the high density res¬ 
onance by Pf-r and adopt the potential, Vo (i?©/r) 3 with 
Vo = 1.25 x 10 -14 eV and the solar radius, R & = 6.96 x 
10 10 cm. Note that Ph is the same for both neutrinos and 
antineutrinos [20] and has an e - sln2 f, i 3 (|Am^ 1 |/£;„) 2/3 v 0 /3 
dependence [21]. Thus, even an order of magnitude un¬ 
certainty in Vo does not have a qualitatively significant 
effect on Ph ■ 

The integrated spectra at SNO and SK: Information on 
the neutrinos emerging from the supernova after the neu- 
tronization burst will be contained in v e and v e spectra 
observed at SNO and SK. 

For the normal hierarchy, the v & flux will be partially 
or completely converted into and v r with the survival 
probability given by [19] 

p = P H P2e + (1 - P H ) sin 2 0 13 . (2) 

The sensitivity of the signal depends on sin 2 20 13 both 
explicitly and implicitly through Ph- The survival prob¬ 
ability for electron antineutrinos is P = Pie [19], which 
is the probability that an antineutrino reaching the earth 
in the V\ mass eigenstate interacts as a 

In the case of the inverted hierarchy, the v e survival 
probability is P = P 2 e [19], which is the probability that a 
neutrino reaching the earth in the V 2 mass eigenstate will 
interact in a detector as v e . The v e survival probability 
is [19] 

P = PnPie + (1 - Ph) sin 2 0 13 . (3) 

Since P\ e and P 2 e depend only on oscillation parameters 
at the solar scale (and the supernova’s zenith angle 0^), 
nothing can be learned about sin 2 20i 3 from the v e (i/ e ) 
flux if the hierarchy is normal (inverted). 

For either hierarchy, we expect little sensitivity to T l , e 
because the survival probability of electron neutrinos is 
no more than about P 2 e ~ sin 2 $12 ~ 0.3. 

For the 32 kton fiducial volume of the reinstrumented 
SK detector (with a 7.5 MeV threshold), and the 1.4 kton 
fiducial volume of the light water tank at SNO (with a 
threshold of 5 MeV), we only consider events that are 
isotropic and indistinguishable from each other; they are 
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FIG. 1. Determination of the binding energy Et,, the super¬ 
nova neutrino mean energies (temperatures) and sin 2 20 i 3 for 
the normal mass hierarchy. The left-hand and right-hand pan¬ 
els correspond to data simulated at sin 2 20 i 3 = 10 -5 (Ph = 1) 
and sin 2 2013 = 10 -2 (Ph = 0), respectively. The cross-hairs 
mark the theoretical inputs, and the 90% and 99% C.L. re¬ 
gions are light and dark shadings, respectively. (E„ x ) is the 
mean energy of the non-electron neutrinos. 

v e + p -» n + e + , (4) 

T O —► P T e — , v e + O —> N + e + . (5) 

We do not consider electron scattering events. The good 
directional capability on these events allows their sep¬ 
aration, and they play an important role in the recon¬ 
struction of the direction of the supernova that in turn 
determines the extent to which earth matter effects may 
be important [22], 

Neutrinos will interact with deuterium in the 1 kton 
fiducial volume of the heavy water tank at SNO (with a 
5 MeV threshold) via the charged current (CC) reactions, 

v e + d —> p + p + e~ , (6) 

v e + d —> n + n + e + . (7) 

In addition we include the reactions of Eq. (5). Two 
neutron captures in addition to a Cherenkov light cone 
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FIG. 2. The same as Fig. 1 but for the inverted hierarchy. 


can distinguish v e -d events from the other charged cur¬ 
rent scattering events on deuterium or oxygen. All NC 
events (whose signal is a single neutron capture and no 
electron), and electron scattering events are neglected. 

To simulate the energy spectra for the channels under 
consideration at the two experiments we assume a typi¬ 
cal supernova [14] at a distance of 10 kpc: E b = 3 x 10 53 
ergs, (E Vc ) : (E De ) : (. E Ux ) :: 0.85 : 1 : 1.1 with 
(E, e ) = 15 MeV, J 7 „ e = 2 (T„ e = (E„ e )/ 3.61), rj*. = 3 
(T Pe = {E P ,)/ 3.99), = 1.5 (T„. = (£„J/3.45), 

L Ve = and L Ve = 1.5 L Vx . We fix sin 2 20 i2 = 0.81, 
Am 2 ! = 5.6 x 10 -5 eV 2 and |Ato3!| = 3 x 10~ 3 eV 2 , 
since variation of these parameters within their future 
bounds has very little effect on the analysis. We gen¬ 
erate v e and v e spectra [22] and then simulate data by 
choosing a point from a Gaussian distribution centered 
at the expectation for the bin and of width equal to its 
square root. In all, there are four spectra; one for SK, 
one for the light water tank at SNO, one for processes 
(5) and (6), and one for process (7). The SK spectrum is 
simulated with 18 bins and the SNO spectra have 13 bins 
each. We simulate four datasets, two for each type of hi¬ 
erarchy and for two values of sin 2 20 i 3 that correspond 
to adiabatic {Ph = 0) and non-adiabatic (Ph = 1) os¬ 


cillation transitions. We perform a y-analysis, freely 
varying E b , (. E Ve ), (E Pe ), ( E „ x ), and sin 2 20 i3 to find the 
90% (Ay 2 < 7.78) and 99% C.L. (Ay 2 < 13.3) allowed 
regions in E b , (E 9tt ), (E Vx ) and sin 2 20i 3 . Although we 
scan in ( E Ve ), we do not count it as a free parameter since 
we do not attempt to determine it (knowing a priori of 
the limited sensitivity to this parameter). We allow the 
ratio L Vt jL Vx to vary between 1 and 2 to accommodate 
both perfect equipartitioning and large departures from 
it. We fix L„ e = L 9b . 

Figure 1 shows the results of this fit for the normal 
hierarchy. The left-hand and right-hand panels corre¬ 
spond to data simulated at sin 2 20 i3 = 10 -5 (for which 
P H = 1) and sin 2 20!3 = 10 -2 (for which Pjj = 0), re¬ 
spectively. In either case, we see that the supernova pa¬ 
rameters can be determined with high precision, but that 
sin 2 2013 is unconstrained. Since the overall normaliza¬ 
tion of the neutrino fluxes depends critically on Eb, it is 
determined with good accuracy. The values of Tp e and 
T Vx are also determined precisely since these parameters 
control the v e spectral distortion (which is independent 
of sin 2 2013) obtained from thousands of events at SK and 
hundreds more at SNO. The experiments are not sensi¬ 
tive to sin 2 2013 because the v e - d events at SNO and u e -0 
events at SK and SNO are statistically insufficient. 

Figure 2 shows the results for the inverted hierarchy. 
Again, the left-hand and right-hand columns correspond 
to data simulated at sin 2 20 i 3 = 10 -5 and sin 2 20i 3 = 
10 —2 , respectively. In the case of non-adiabatic transi¬ 
tions, an upper bound on sin 2 20i 3 can be placed. For 
adiabatic transitions, sin 2 20 i 3 and (E^ e ) cannot be si¬ 
multaneously bounded if both are left free. When we 
restrict {Ep e ) to lie between 10.5 MeV and 19 MeV 
and (E Vx ) / (Ea e ) to be larger than about 0.7 (recall that 
(E„ x ) / {Ev e ) is expected to be larger than unity), a lower 
bound on sin 2 20i 3 is obtained. For values of sin 2 20 i 3 
between 10~ 4 and 10~ 3 , an upper or lower bound can 
be placed on sin 2 20 i 3 depending on whether it is closer 
to 10” 4 or to 10^ 3 . In the inverted hierarchy, we have 
a lesser sensitivity to Tp e and T Vx since the v e flux is 
also sensitive to sin 2 20 i 3 leading to competition between 
these parameters. Bounds on sin 2 20 i 3 can be placed be¬ 
cause the u e spectrum is more sensitive to sin 2 20 i 3 in the 
inverted hierarchy and the v e signal at SK is huge. 

Although the regions in Figs. 1 and 2 are calculated 
assuming that the neutrinos detected at SK crossed both 
the mantle and core of the earth (cos0z = —0.93), and 
those at SNO crossed the mantle only (cos0^ = —0.1), 
we have established that the bounds placed are largely 
independent of the supernova’s zenith angles at the two 
experiments. 

If the mass hierarchy is unknown at the time of a su¬ 
pernova signal, it can be deduced provided sin 2 20 i 3 > 
10~ 3 [19,23], and the hierarchy is inverted. For values of 
sin 2 2013 smaller than « 10 -4 , Ph is not close to zero and 
the survival probabilities are similar for the two hierar- 
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dries rendering them indistinguishable [19]. In the case 
of a normal hierarchy, we see from Fig. 1 that sin 2 29\ 3 is 
unconstrained even for adiabatic transitions, thereby in¬ 
dicating a lack of discriminatory power between Ph = 0 
and Ph = 1 or equivalently between the mass hierar¬ 
chies. On the other hand, for the inverted hierarchy and 
adiabatic transitions, a lower bound on sin 2 29 33 can be 
placed which in turn means that the inverted hierarchy 
can be selected over the normal hierarchy. 

Future prospects: The next generation of proton decay 
experiments such as Hyper-Kamiokande (HK) [24] and 
UNO [25] are expected to offer a new level of sensitivity 
to the physics of supernovae and neutrino mixing. We 
consider the proposed 1 Mton HK detector. With no 
specific information about the detector, we treat it as a 
scaled-up version of SK. We assume a fiducial volume for 
supernova neutrinos of 890 kton, which is consistent with 
the fiducial volume to total volume ratio expected for the 
proposed UNO detector [26]. 

We find that our qualitative conclusions for SK and 
SNO continue to hold for HK. The quantitative differ¬ 
ences are easily anticipated as a result of its larger vol¬ 
ume. The supernova parameters can be determined with 
greater accuracy although T„ e will remain unknown. In 
the case of the inverted hierarchy and adiabatic transi¬ 
tions, while T„ x can be determined without theoretical 
prejudice, a plausible window has to be chosen for T Pe to 
constrain sin 2 29 33 . Also, tighter upper or lower bounds 
can be placed on sin 2 26 * 13 . We emphasize that in the case 
of a normal hierarchy, both sin 2 29 13 and the hierarchy 
remain unknown. 

Summary: We have considered what information can be 
extracted from neutrinos detected at SNO and SK from 
a galactic supernova. The information they carry is of 
major importance in understanding the astrophysics of 
supernovae. The binding energy released in the super¬ 
nova and the temperatures of the non-electron neutrinos 
expelled may be determined with good precision for most 
values of sin 2 26 * 13 . Bounds on sin 2 29 13 can be placed if 
the neutrino mass hierarchy is inverted. In this case the 
hierarchy can be determined if sin 2 29 13 > 10 -3 . 

The above conclusions apply to Hyper-Kamiokande as 
well. 
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